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ABSTRACT 

Molecular clouds (MCs) are highly structured and "turbulent". Colliding gas streams of atomic 
hydrogen have been suggested as a possible source of MCs, imprinting the filamentary structure as 
a consequence of dynamical and thermal instabilities. We present a 2D numerical analysis of MC 
formation via converging HI flows. Even with modest flow speeds and completely uniform inflows, 
non-linear density perturbations as possible precursors of MCs arise. Thus, we suggest that MCs 
are inevitably formed with substantial structure, e.g., strong density and velocity fluctuations, which 
provide the initial conditions for subsequent gravitational collapse and star formation in a variety of 
galactic and extragalactic environments. 

Subject headings: turbulence — methodsmumerical — ISMxlouds — ISM:kinematics and dynamics — 
stars:formation 
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1. MOTIVATION 

Molecular clouds (MCs) in our Galaxy are com- 
plex and highly-structured, with broad, non-thermal 
linewidths suggesting substantial turbulent motion 
(jFalgarone &_Phihps 1990; WiUiams et al. 2000). Thus, 
MCs very likely are not static entities and might not 
necessarily be in an equilibrium state, but their proper- 
ties could well be determined by their formation process. 
The importance of initial conditions for cloud structure is 
emphasized by observational and theoretical evidence for 
short cloud " lifetimes" ( Ballestcros-Parcdes ct al. 1999; 
lElmegreenl200nHHartmann et al.-200L:„Hartmann.,2002L 

Flows are ubiquitous in the interstellar medium (ISM) 
due to the energy input by stars - photoionization, winds, 
and supernovae. In principle, they can pile up atomic gas 
to form MCs. Shock waves propagating into the warm 
ISM will fragment in the p resence of thermal instabil- 
ity and linear perturbations llHennebelle fc Peraumi999l 
IKova,ma, Tnutsukal [20001 1200?. '2004b^ a,nd a,1- 
low H2-formation within a few Myrs ( Bcrgin ct al. 200^ 
in a plane-parallel geometry. We envisage the collid- 
ing flows less as e.g. interacting shells, but as (more 
or less) coherent gas streams from turbulent motions on 
scales of the order of the Galactic gaseous disk height 
((Ballesteros-Parcd es et al.llT9 99: Hartmann ct al. 2001). 
Parametrizing the inflows as a ram pressure allowed 
iHennebelle et all l|2003l l2004|) to study the fragmenta- 
tion and collapse of an externally pressurized slab. 

In this paper we present a study of the generation of fil- 
aments and turbulence in atomic clouds ~ which may be 
precursors of MCs - as a consequence of their formation 
process, extending the model of large-scale c olliding HI- 
flows outlined by IBallesteros-Paredes et alJ |.1999) and 
IBurkerHll200l . We discuss the dominant dynamical and 
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thermal instabilities leading to turbulent flows and frag- 
mentation of an initially completely uniform flow. Re- 
sulting non-thermal linewidths in the cold gas phases 
(the progenitors of MCs) reproduce observed values, em- 
phasizing the ease with which turbulent and filamentary 
structures can be formed in the ISM. 

This study is a " proof of concept" , discussing the 
dominant instabilities leading to non-linear density per- 
turbations and quantifying the timescales that are nec- 
essary to reach conditions under which molecules - 
and eventually stars - could form. Since turbulence 
is thought to be one of th e main agents controlling 
star formation (Larson 198lt iMac Low & Klessenll2004l: 
lElmegreen fc Scala2004) , and consequent density fluctu- 
ations surely are crucial to gravitational fragmentation, 
star formation theory can benefit from a better under- 
standing of the structure initially present in MCs. De- 
tailed aspects and pa rameter studies will be presented in 
a forthcoming paper l)IIeitsch et al.ll2005j) . 

2. PHYSICAL MECHANISMS 

We restricted the models to hydrodynamics including 
thermal instabilities, leaving out the effects of gravity 
and magnetic fields. Gravity would eventually lead to 
further fragmentation, and magnetic fields would be ex- 
pected to have a stabilizing effect. For this regime, then, 
we identify three relevant instabilities: 

(1) Th e Non-linear Thin Shell Instability (NTSI, 

arises in a shock-bounded slab, when rip- 
ples in a two-dimensional slab focus incoming shocked 
material and produce density fluctuations. The growth 
rate is ~ Csfc(fcA)^/^, where Cs is the sound speed, k is 
the wave number along the slab, and A is the amplitude 
of the spatial perturbation of the slab. Numerical stud- 
ies focused on the generation of substructure via Kelvin- 
Helmholtz- modes jB londin & M ark£lll99i|) . on the role 
of gravity llHunter et al. nl986|) and on the effect of the 
cooliM stre ngth l|Hueckstaedd l200l . iWalder fc Folinil 
11998', '2000*) discussed the interaction of stellar winds, 
and Klein & Woods ( 1998) investigated cloud collisions. 
The latter authors included magnetic fields, albeit only 
partially. 

(2) The flows deflected at the slab will cause Kelvin- 
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Helmholtz Instabilities (KHI), which have been studied 
at great length analytically and numerically. If the veloc- 
ity profile across the shear layer is given by a step func- 
tion, and if the densities are constant across the layer, 
the growth rate is kAU, where AU is the velocity dif- 
ference. If aligned with the flow, magnetic tension forces 
can stabilize against the KHI. 

(3) The Thermal Instability fTI. lFieldllT965|) rests on 
the balancing of heating and cooling processes in the 
ISM. The TI develops an isobaric condensation mode 
and an acoustic mode, which - under ISM-conditions 
- is mostly damped. The condensation mode's growth 
rate is independent of the wave length, however, since 
it is an isobaric mode, sm aller perturbations will grow 
first ijBurkert fc LinI 1200 (^1. A lower growth scale is set 
by heat conduc tion, whose scale needs to be resolved 
ijKovama fc Inu tsuka 2004a). The signature of the TI are 
fragmentation and clumping as long as the sound cross- 
ing time is smaller than the cooling time scale. The TI 
can drive tvirbulcncc in an otherwise quiescent medium 
ijAudit fc H cnncbcUc 2005) - even continuously, if an 
episodic heating source is available ijKritsuk fc NormanI 

3. NUMERICAL METHOD 

All three instabilities grow fastest or at least first 
on the smallest scales. This poses a dire challenge for 
the numerical method. We c hose a metho d based on 
the 2nd order BGK formalism i j Prend crgast 
i Slvz fc PrendergastllTOQi iHeitsch et al.ll2004fls!^reral] 
l2005j) ■ allowing control of viscosity and heat conduction. 
The statistical properties of the models are invariant un- 
der changes of viscosity, heat conduction and grid reso- 
lution, although the flow patterns change in detail — as 
to be expected in a turbulent environment ('Hcitsc h et al.l 
|2005). The linear resolution varies between 512 and 2048 
cells. The heating and cooling rate s are restricted to op - 
tically thin atomic lines following iWolfire et al.l l) 19951) . 
so that we are able to study the precursors of MCs 
up to the point when they could form H2. Dust ex- 
tinction becomes important above column densities of 
A^(III) « 1.2 X lO^^cm"^, which are reached only in 
the dense cold regions of the models. Thus, we can use 
the unattenuated U V radiation field for grain heating 
(|Wolfire et alJlT995l) for most of the simulation domain, 
expecting substantial uncertainties in cooling rates only 
for the densest regions. The ionization degree is derived 
from a balance between ionization by cosmic rays and 
recombination, assuming that Ly a photons are directly 
reabsorbed. 

Two opposing, uniform, identical flows in the x-y com- 
putational plane initially collide head-on at a sinusoidal 
interface with wave number ky = 1 and amplitude A. 
The incoming flows are in thermal equilibrium. The sys- 
tem is thermally unstable for densities 1 ^ n < lOcm"'^. 
The cooling curve covers a density range of 10~^ ^ n < 
10^ cm~'^ and a temperature range of 30 < T < 1.8 x lO'' 
K. The box side length is 44pc. Thus, Coriolis forces 
from Galactic rotation are negligible. For a interface with 
A = 0, a cold high-density slab devoid of inner structure 
would form. The onset of the dynamical instabilities thus 
can solely be controlled by varying the amplitude of the 
interface perturbation. This allows us to test turbulence 
generations under minimally favourable conditions. 



Model series B has uq — 1.0cm and Tq = 2.5 x 
lO^K, series C has no = 0.5cm~3 and Tq = 5.3 x lO^K. 
The amplitude of the interface perturbation corresponds 
to 2.5% of the box length, i.e. l.lpc in the horizontal 
direction. The second digit of the model name gives the 
Mach number of the inflow. The inflow velocity varies 
between 6 < Wi„ < 17 km s^^. 

4. RESULTS 

We chose to present a few models containing the salient 
features of cloud formation via colliding flows. We ask 
how hard it is to generate non-linear turbulent density 
perturbations in an otherwise uniform flow. 

4.1. Dominant Instabilities 

The structures generated in colliding flows depend 
strongly on the initial parameters (Fig.^^, as a result of 
the dominating instability. This is not surprising, since 
all three instabilities at work have different signatures. 
For high-density, low- velocity inflows (model Bl, upper 
left panel), the TI dominates and leads to fast cooling, 
manifested in a coherent slab of cold gas. This situation 
comes closest to the ID plane-parallel slab. 




Fig. 1. — Temperature maps (logT in K) for models represen- 
tative of the dominating instability (see text). From upper left to 
lower right: Bl (TI), CI (KHI), C2, (KHI+NTSI), B3 (NTSI). 

Reducing the density (model CI, upper right) leads 
to slightly less efficient cooling, thus giving the dy- 
namical instabilities time to work, in this case domi- 
nated by the KHI. The eddies are visible at the flanks 
of the initial slab. Increasing the inflow speed while 
keeping all other parameters fixed increases the verti- 
cal x-momentum transport, so that for model C2 (lower 
left), the NTSI will arise. Its typical signature are long 
strands of denser (and colder) gas predominantly along 

All plots without time-dependence are taken at the endpoint 
of the corresponding model, i.e. at ISMyrs for model B3, and at 
19Myrs for all other models. 
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the flow direction (Hueckstaedt 2003). The KHI modes 
that are still discernible here have vanished in model 
B3. Here, the NTSI dominates the dynamics almost 
completely. Oblique colliding flows (not shown) lead 
to a nearly instantaneous break-up of the initial slab, 
because they excite KHI modes at the scale of the ini- 
tial, non-linear perturbation. The NTSI dominates more 
and more with higher Mach numbers (up to Mach 4, not 
sh own). This morpholo gical discussion will be quantified 
in lHeitsch erall l|2005|) . 

Despite the symmetric initial conditions, all models de- 
velop large-scale asymmetries, resulting from an amplifi- 
cation of slight differences at machine accuracy between 
the upper and lower half of the domain. The main culprit 
is the strong cooling. Without cooling (i.e. for an adia- 
batic or isothermal equation of state), or for ID problems 
with cooling, the code preserves perfect symmetry. 

4.2. Mass Distribution 

Figure |21 shows the mass spectrum of cores. Cores 
are deflned as coherent regions with densities n > 
lOOcm^^. Each of the histograms comprises approxi- 
mately 50 cores. The dashed line denotes a spectral in- 
dex of —1.7, as observe d for observed molecular cores 
ijHeithaus en et af] ll99^. The resolution limit for both 
histograms lies at 1.8Mq. Note that (a) these "cores" 
correspond to cold HI regions, not molecular cores, and 
that (b) the masses are per length in our 2D models. 
A direct comparison between these 2D spectra and ob- 
served 3D mass spectra requires assumptions concern- 
ing what structures the fllaments would correspond to in 
three dimensions. Assuming that the probability density 
functions are the same for 2D and 3D, the mass spectr a 
are expected to be flatter in 3D ijChappell fc Scalol200H) . 
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Fig. 2. — Left: Mass-histogram of objects with densities n > 
lOOcm"^, i.e. likely precursors of molecular fragments, for models 
C2 and B3. The dashed line denotes a power law with exponent 
— 1.7. The resolution limit is Mr = I.SMq. Right: Histogram of 
the local Jeans number L/Xj for the cold regions, for models C2 
and B3. Clearly, none of the dense objects would be gravitationally 
bound. 

Would these cold cores be gravitationally unstable? 
Although the models do not include self-gravity, we can 
estimate the thermal Jeans length Aj = (7r/(Gp))^/^Cs 
and its turbulent counterpart \t = {n/ {Gp)y/^ {v^y^^ . 
The ratio of the cores size over the core's Jeans length 
nj = L/Xj < 1, i.e., none of the cores would be 
gravitationally unstable (Fig|21 right). Since \j oc 
T^l"^ jiPjTY^'^ oc T for an isobaric contraction, the his- 
togram of Figure El (right) would shift by a factor of 
3 to larger Jeans masses if we cooled the gas down to 
T w lOK, thus still yielding Jeans-stable objects. The 
core size is determined by the geometric mean of the 
longest and shortest radius. 

However, if we determine the "global" Jeans number of 
the cold gas — given by the thickness of the slab over the 
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Fig. 3. — Average length scale (thick lines), turbulent Jeans 
length (medium lines) and thermal Jeans length (thin lines) for 
models C2 and B3 against time. 



global Jeans length derived from the mean density and 
temperature in the cold gas — after lOMyrs for model 
B3 and 15Myrs for model C2, the length scale of the cold 
gas L > Xj, i.e. the "slab" would become gravitationally 
unstable without turbulence (Fig. |3J). Note that these 
quantities are global measures in the sense that they do 
not refer to isolated cold regions. The turbulent Jeans 
length Xt > L for all times and models. Global (i.e. 
large-scale) gravitatio nal effects c ould still have a crucial 
effect on the system llBurkert fc Hartman'nll200^ . espe- 
cially once the inflow stops. 

4.3. Linewidths and Kinetic Energy Modes 

A primary observable of interstellar clouds is the line- 
of-sight velocity dispersion a-u. MCs co nsistently show 
non-t hermal linewidths of a few km s~^ l|Williams et alJ 
12000(1 ■ that - together with temperatures of T ?s lOK 
- are generally interpreted as supersonic turbulence in 
those clouds. The linewidths in our models are consis- 
tent with the observed values (Fig. Q . The " observed" 
linewidth is derived from the density- weighted histogram 
of the line-of-sight velocity dispersion in the cold gas at 
T < lOOK (flUed symbols). Comparing this to the inter- 
nal linewidth of coherent cold regions (open symbols), we 
note a marked offset between the two values. The sound 
speed of the cold gas ranges around 0.8km s~^. Thus, the 
internal velocity dispersions do not reach Mach numbers 
> 1. Hence, the "supersonic" linewidths are a conse- 
quence of cold regions moving with respect to each other, 
but not a result of internal supersonic turbulence in the 
cold gas which eventually would be hosting star forma- 
tion. T his is consistent with the argument bv lHartmannI 
l|2002l) . that because of the ages and small spatial disper- 
sions of young stars in Taurus, their velocity dispersions 
relative to their natal gas are very likely subsonic. Note 
that the turbulent linewidths amount only to a frac- 
tion of the inflow velocity. 

Figure El shows the compressible, the solenoidal and 
the total specific kinetic energy for the whole domain 
(left) and for the cold gas {T < lOOK, right). Because of 
the highly compressible initial conditions, the total spe- 
cific kinetic energy is dominated by compressible modes. 
However, within the bounding shocks, the highly com- 
pressible inflows are converted efficiently into solenoidal 
motions. Note that because of the 2D geometry, the ra- 
tio of solenoidal over compressible kinetic energy is only 
a lower limit. An extension of Figure|Slto radiative losses 
and internal energy allows us to estimate the overall ef- 
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Fig. 4. — Mean velocity dispersion in cold gas at T < lOOK (filled 
symbols) and in the cold coherent regions (open symbols) against 
the inflow velocity. 



ficiency of turbulence generation in MCs ijHeitsch et al.l 
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Fig. 5. — Specific kinetic energy fraction against time for models 
C2 and B3, split up into compressional, solenoidal and total specific 
kinetic energy. Left: For the whole domain. Right: for the cold 
gas (T < lOOK). The specific kinetic energy is normalized to the 
value at t = 0. 



5. SUMMARY 

Even for completely uniform inflows, we have shown 
that the combination of dynamical and thermal insta- 
bilities efhciently generates non-linear density perturba- 
tions that seed structure of eventual MCs. There is a 



direct correlation between the morphology of the result- 
ing clouds and the dominating instability. 

While our "cloud" would be gravitationally unstable, 
the isolated cold regions would still be stable against 
gravity. Fragmentation and turbulent mixing driven by 
the incoming warm gas would prevent the global col- 
lapse of the cloud - under the caveat that a finite extent 
of the cloud in the vertical direction might lead to edge 
effects resulting in collapse (Burkcrt & Hartmann 2004). 
Linewidths reached in the cold gas are consistent with 
observed values of a few km s~^ (Fig 01 and reach only 
a fraction of the inflow speed. The internal linewidths, 
however, are generally subsonic. Thus, the label "super- 
sonic turbulence" refers to the velocities with respect to 
the cold gas, but does not necessarily give a hydrody- 
namically accurate description of the cold gas. 

Although we adopted a specific cooling curve and thus 
set the physical regime for our models, we expect the 
mechanism to work on a variety of scales. The surface 
density of the cold gas should give us a rough estimate of 
the amount of stars forming later on. Even though the 
cold gas mass depends strongly on the turbulent evolu- 
tion of the slab, it correlates strongly with the inflow 
momentum. In this sense, colliding flows not only could 
explain the rather quiescent star formation events as in 
Taurus, but would be a suitable model for generating 
star bursts in galaxy mergers. 
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